Mean-field phase diagram of the 1-D Bose gas in a disorder potential 
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We study the quantum phase transition of the ID weakly interacting Bose gas in the presence of 
disorder. We characterize the phase transition as a function of disorder and interaction strengths, 
by inspecting the long-range behavior of the one-body density matrix as well as the drop in the 
superfiuid fraction. We focus on the properties of the low-energy Bogoliubov excitations that drive 
the phase transition, and find that the transition to the insulator state is marked by a diverging 
density of states and a localization length that diverges as a power- law with power 1. We draw 
the phase diagram and we observe that the boundary between the superfiuid and the Bose glass 
phase is characterized by two different algebraic relations. These can be explained analytically by 
considering the limiting cases of zero and infinite disorder correlation length. 
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I. INTRODUCTION 

The effect of disorder on quantum systems is a subject 
of both fundamental and practical interest. Since the 
seminal work of Anderson it has become clear that 
what seems at first sight a nuisance is in fact a source 
of very rich physical behavior: the disorder is not just a 
perturbation of the wave functions, but localizes the low 
energy states in three dimensions; in lower dimensions 
even all single particle states are localized. 

A single particle picture can be a good approxima- 
tion of a dilute gas, but in practice interactions be- 
tween particles play often a crucial role. Understand- 
ing of the interplay of disorder and interactions is there- 
fore of fundamental importance. In the context of solid 
state physics, the fermionic problem is the most relevant 
one, but thanks to the enormous progress in experimen- 
tal control over ultracold atomic gases, the problem of 
the disordered Bose gas has become a subject of a vigor- 
ous research activity as well. Anderson localization has 
been experimentally observed recently in one dimensional 
bosonic systems with vanishing interaction [2|, |3| • More 
recent experiments have moved away from the limit of 
vanishing interactions, so to study the localized (Bose 
glass) to superfiuid phase transition 0, [1[ . 

The theoretical interest in this phase transition dates 
back much longer and a variety of theoretical techniques 
have been used to tackle the problem. There are two 
main regimes that have been considered, one marked by 
weak disorder and arbitrary interactions, the other char- 
acterized by weak interactions and arbitrary disorder. 
The former has been the object of the first investigations 
in the eighties. Using a renormalization group analysis, 
Giamarchi and Schultz Q were able to study the quan- 
tum phase transition in the limit of weak disorder in one 
dimension. The picture emerging from this analysis was 
that for a finite amount of disorder, a minimal strength 
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of interactions is required to break the Anderson localiza- 
tion, but that for too strong interactions, the system is 
driven into a strongly correlated localized Bose gas phase. 
Their renormalization group approach was able to study 
the latter phase transition quantitatively, but it could not 
be clarified whether the transition at the weak interaction 
side is of the same nature. The interplay between peri- 
odic and disordered potentials was first addressed in the 
seminal work by Fisher et al. [7| , where the insulating dis- 
ordered phase, named 'Bose glass', was contrasted to the 
Mott insulator phase by its compressible nature and to 
the superfiuid phase by its vanishing superfiuid stiffness. 
Quantum Monte Carlo Q and density matrix renormal- 
ization group 0] studies have investigated in detail the 
disorder phase diagram in the limit of strong interactions 
close to the Mott insulator phase. 

Early experimental efforts to reach the Bose glass 
phase coming from the Mott insulating phase were pre- 
sented in Ref. 10]. More recently several experimental 
works have addressed the quite different regime where 
the gas is in the weakly interacting limit and with many 
particles in each potential minimum. Above a critical 
strength of the disorder, fragmentation of the conden- 
sate density and loss of spatial coherence were identi- 
fied 0, IE)]- All these experiments have considered the 
one dimensional geometry. We will restrict our analy- 
sis to this case as well. Thanks to the weak interaction 
limit, the theoretical description of these experiments can 
be performed in a first approximation with the Gross- 
Pitaevskii equation. The Gross-Pitaevskii equation in 
the presence of a disorder potential was recently studied 
by several groups [ll|, [H, [H, [lij]- For the condensate 
wave function, two different cases were identified: a con- 
nected density profile at weak and a fragmented one at 
strong disorder. The properties of the elementary Bo- 
goliubov excitations stemming from the Gross-Pitaevskii 
state were studied analytically in the superfiuid regime. 
Their localization length was shown to exhibit a power 
law behavior as a function of energy: E a , with a — 2 
deep in the superfiuid phase 13| and a — 1 at the phase 
transition [14j |. Numerically we found that a < 1 in 
the Bose glass phase [15[. The density of states (DOS) 
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was shown to be constant in the superfluid phase as for 
phonons in random chains [16j |: in the Bose glass phase, a 
divergence of the low energy DOS was numerically iden- 
tified in our previous work. This result is in accordance 
with the real space renormalization group analysis by 
Altman et al. 17] with particle hole symmetry. Fisher 
et al. Q on the other hand argued that the low energy 
DOS should be constant in the Bose glass phase. This 
difference in behavior could be due to a different nature 
of the phase transition between the superfluid and glassy 
phase for weak and strong interaction as suggested by 
Giamarchi and Schultz Q. 

In the present work, we make a deeper analysis of the 
condensate wave function and Bogoliubov excitations in 
the different phases. We carry out a new analysis of 
the superfluid fraction across the phase transition. In 
Ref. |l5|, the phase diagram was restricted to the limit 
of weak interaction energy where the healing length £ 
was large as compared to the correlation length of the 
disorder r\. In that case, the disorder is effectively an 
uncorrelated white noise (WN) potential. In this work, 
we extend the phase diagram to the regime where the 
r\ 3> £. Kinetic energy is then not important to determine 
the density profile of the condensate and the Thomas- 
Fermi (TF) approximation is accurate. The WN and TF 
regimes are marked by two different power-law relations 
between interaction and disorder at the phase boundary. 
For these we provide a rigorous analytical arguments. 

The paper is organized as follows. In Sec. [II] the the- 
oretical model is presented in detail. Sec. IIIII is devoted 
to the predictions of the extended mean-field model. In 
Sec. IIVI the study of the correlation length is presented. 
An analysis of the density of states and localization prop- 
erties of the Bogoliubov excitations is carried out in Sec. 
IVl In Sec. I VII a study of the superfluid fraction of the 
gas is performed. The phase diagram is discussed in Sec. 
IVIII Our conclusions are presented in Sec. I Villi 



II. THEORY 

We want to study the properties of the ID Bose gas at 
zero temperature. These gases are well described by the 
mean field theory in 3D, but in lower dimensions no real 
condensate is present and the standard mean field theory 
is no longer valid. Nevetheless, the Bogoliubov approach 
can be extended to weakly interacting low-dimensional 
bosonic systems in a density-phase version defined on a 
lattice [18[ ■ A homogeneous ID Bose gas is in the weakly 
interacting regime when p£ >• 1, where £ = h/^/mpg 
is the healing length, p the total density, g the coupling 
constant and m the mass. 

The Hamiltonian describing the bosonic system is 



H = dr 



,(1) 



where Hq — — h 2 d^/(2m) + V(r) is the single-particle 
Hamiltonian, i& is the field operator and V(r) the exter- 



nal potential. Here we study the case where V(r) is a 
Gauss-distributed and Gauss-correlated disorder: 



(V(r)V(r')) 



A 2 e 



(r—T 1 ) 



(2) 



where A is the disorder amplitude and rj is the spatial 
correlation length. We make the general choice of Gauss- 
distributed disorder. We are aware that many experi- 
ments aimed at the characterization of the phase tran- 
sition deal with speckle potentials @, These poten- 
tials have a lower bound and their statistical distribution 
does not show a lower Gaussian tail. Our analysis could 
be easily extended to this V(r) is only treated 

numerically. On the other hand, the assumption of spa- 
tial Gauss-correlation is consistent with the experimental 
realizations [5|. 

In low dimensionality, the Bogoliubov approach re- 
quires the definition of the field in terms of density p 
and phase 9 operators as = e l6 \fp. In the mean field 
approach one splits the density operator in a c-field part, 
Po, and a fluctuation term 8 p. A correct definition of the 
phase operator is only possible fl8j by introducing a spa- 
tial discretization of step size £. Provided that the system 
is in the high density regime, pi = (po + (Sp))t > 1, and 
that every cell is largely populated, the phase operator 
can be defined precisely. The density fluctuation can then 
be treated perturbatively: this perturbative approach is 
valid under the assumption of small density fluctuations, 
Sp/p <C 1, and spatially slowly varying phase fluctua- 
tions, 9 i+ i 

The ground state density profile po obeys the Gross- 
Pitaevskii equation 



Ho+gpo{r) y/ po{r) = py/p {r), 



(3) 



where p is the chemical potential. The excitations 
of the system can be computed via the Bogoliubov-dc 
Gennes equations, obtained by linear expansion around 
the Gross-Pitaevskii solution 



H +2gp (r)- pjuj(r)+gp (r)vj(r) = EjUj(r), 

9Po( r ) u ]( r ) - [Ho + 2gPo(r) - pj vj(r) = E 3 Vj(r), 

(4) 



that define the Bogoliubov uj(r) and Vj(r) modes, nor- 
malized as 



drK-(r)| a -|« 3 -(r)| 2 = l. 



(5) 



Phase and density operators can be expressed in terms 
of these wave functions as 11911 



1 



Sp(r) = Jpjr)^^ + p*{r)P], 



(6) 
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where bj and bj are bosonic annihilation and creator op- 
erators of an excitation with energy Ej and with 

P.i(r) + 0j(r) Ps(r)-6j(r) 



(7) 



The extended Bogoliubov approach requires an orthog- 
onalization of these modes with respect to the ground 
state [HI], which finally brings to the new modes u±j{r) 
and v±j(r). 

Taking advantage of Wick's theorem, the one-body 
density matrix at T = 0, computed with the extended 
Bogoliubov method, takes the form 



G(r,r') = ^p(r)p(r>)t 



3 



v±j(r) 



(r') 



(8) 

where the only contribution comes from the quan- 
tum fluctuations. This two-point function is not self- 
averaging. We will preferably work with the spatially 
averaged degree of coherence 



9i(r) = \ 



dr'- 



G(r,r') 
y/p(r)p{r') 



(9) 



that directly gives information about the decay of the 
one-body density matrix. This quantity, in the quasi- 
condensed phase, is expected to show an algebraic de- 
cay 0, as in the homogeneous case [loj], together with 
a finite superfluid fraction. In presence of disorder the 
system can undergo a quantum phase transition to the 
Bose glass phase characterized by an exponential decay 
of the one-body density matrix 



III. GROUND STATE AND PHASE 
FLUCTUATIONS 
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FIG. 1: (color online) Quasicondensed phase: a. ground state 
wave function and b. first two excitations v±j for U — 1.44 E c 
and A = 0.8 E c . 
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Three energy scales enter our problem: the amplitude 
of the disorder A, the interaction energy, U = gNo/L, 
where N is the number of particle in the ground state, 
and the energy associated to the correlation length of 
the potential E c — 2 . In the analysis that follows 
the problem is rescaled with respect to this latter. We 
have performed numerical calculations on systems of fi- 
nite size (up to 4096^), adopting periodic boundary con- 
ditions for the equations and for the disordered potential 
V(r). Every quantity has been averaged on several disor- 
der configurations. By varying the system size, we could 
extract the limiting behavior of each quantity in the ther- 
modynamic limit. In order to describe the system in the 
continuous limit wc fulfill the condition that the kinetic 
hopping energy t = be much larger than any other 
characteristic energy of the system. This must hold in 
particular for the energies U , A, p and E c . In the numer- 
ical simulations that follow we took an average quasicon- 
densate density of Norj/L = 8. It is important to remark 



FIG. 2: (color online) Bose glass phase: a. ground state wave 
function and b. first two excitations v±j for U = 0.48 E c and 
A = 0.8 i? c . The thin grey line represents an excitation at 
higher energy, namely v±4,o(r). 



that in the mean field limit the phase boundary, deter- 
mined from the Bogoliubov-de Gennes equations and the 
functional shape of the degree of coherence, does not de- 
pend on the density and interaction strength separately, 
but only on the product of the two quantities that ap- 
pears in Eqns. and ([!]). As outlined in Sec. [Til this 
model is valid in the limit of large density and its predic- 
tion are increasingly accurate in the limit p — > 00, g — > 
at constant gp. 

In Fig. [T]a the ground state is shown for A = 0.8 E c 
and U = 1AAE C1 in the quasicondensed phase, computed 
via Eq. ([3]). In Fig. [T]b the first two excitations v±j(r) 
appear (J = 1,2). The same quantities computed in the 
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Bose glass phase (A = 0.8 E c and U = 0.48 E c ) are shown 
in Fig. H 

It is evident, comparing Figs. [T]a and[2la, that in the 
case of stronger interaction the ground state wave func- 
tion is smoother. In fact the interaction term involves a 
more homogeneous effective potential in Eq. (|3|) where 
the bare disorder potential is screened. In a case dom- 
inated by the disorder, the ground state separates into 
fragments linked by regions with an exponentially van- 
ishing wave function, that can be seen as weak links. 

Low energy excitations have a delocalized shape and 
follow the modulation of the ground state wave function, 
but, differently from the ground state, they show a phase 
character with a number of nodes that increases for in- 
creasing energy. In the homogeneous case the lowest en- 
ergy excitations are found to be plane waves. The phase 
fluctuations preserve a plane wave profile in the quasi- 
condensed phase, as can be seen in Fig. [TJb, only slightly 
modulated by the underlying disordered potential. On 
the other hand, at smaller values of U, in the Bose glass 
case, the disorder starts to compete with the interaction 
and the v±j(r) modes start losing their regular shape, de- 
veloping nodes in correspondence with low density zones, 
i.e. with high barriers of the potential. This is shown in 
Fig. Ob, together with an example of an excitation at 
higher energy that displays a fastly oscillating behavior 
and does not contribute in determining long-range prop- 
erties. 

It is instructive to analyze the role of the correlation 
length rj on the ground state and on the excitations. The 
quantity E c and its relation with // turn out to be essen- 
tial in determining the regime of the Bose gas. In fact 
small values of U/E c and A/E c imply that r\ <C £ and 
the ground state wavefunction is spread over many corre- 
lation lengths of the potential. In this limit, the disorder 
is equivalent to a WN potential. On the contrary, for 
large values of U/E c and A/E c , the system enters the 
TF regime, where the kinetic term is negligible and the 
ground state follows the spatial variations of the poten- 
tial. In Fig. [3^ the ground state is shown for a case 
close to the WN limit (U = A = 1.6 X lO" 3 ^). The 
same quantity is shown in Fig. 0^ for U = A = 25.6 E c , 
close to the TF regime. The differences explained above 
are evident in the distributions of the ground state that 
is spread over many correlation lengths in the WN case, 
whereas it fills the potential minima in the TF regime 
and it approaches the form 



for 



2 _ a - Vjr) 

UL ' 
1 — 11 otherwise. 



\Mr)\ 



\Mr)V = 0, 



fi > V{r) 



(10) 



The WN case presented in Fig. [3] is in the superfluid 
phase. The low lying excitations are very close to plane 
waves with small modulations and a regular spacing be- 
tween the nodes. The TF example (see Fig. [4} is in 
the Bose glass phase and the nodes of the excitations 
are pinned to the low condensate density regions. Their 
profile follows closely the condensate amplitude. 
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FIG. 3: (color online) a. Ground state and b. first two 
excitations v±j for U = A = 1.6 x 10~ 3 E c , where the system 
is in the superfluid phase and in the WN regime. 
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FIG. 4: (color online) a. Ground state and b. first two 
excitations v±j for U = A = 25.6 E c , where the system is in 
the Bose glass phase and in the TF regime. 



Note that the ratio between interaction and disorder 
amplitudes is in both cases A/U = 1. The transition 
between superfluid and Bose glass phases can apparently 
be tuned by only varying the disorder correlation length 
77. We will come back to this point in Sec. I VIII where we 
present the phase diagram. 



IV. CORRELATION LENGTH 

As outlined above, the phase transition is characterized 
by inspecting the long range behavior of the degree of 
coherence expressed by Eq. ©. Fig. [5] shows the degree 
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FIG. 5: (color online) Degree of coherence gi(r) for fixed value 
of disorder A = 0.8 E c for different values of the interaction 
energy. The cases with U = 0.48 — 0.8 E c are in the Bose glass 
phase, whereas for U = 1.12 — 1.92 E c are superfluid. 




100 200 300 400 500 

Distance [ri] 

FIG. 6: (color online) a. G(ro, r) in the superfluid and in the 
Bose glass phases, b. corresponding low-energy v± excita- 
tions in the Bose glass phase. The arrows point out the link 
between the jumps in G(ro, r) and the nodes of the low-energy 
excitations. 



of coherence for the same value of disorder (A = 0.8 E c ) 
varying the interaction energy. We show only the spatial 
interval [0, L/A] that reflects the L — > oo behavior: at 
longer distance, deviations due to the periodic boundary 
conditions affect this quantity, as verified by studying the 
scaling with L. For the lowest value of the interaction 
energy, g\ (r) shows an exponential decay that marks the 
Bose glass phase in a situation where the disorder breaks 
the long-range coherence. By increasing the interaction 
energy the decay becomes slower up to a certain value of 
interaction that drives the system in the quasicondensed 
phase characterized by an algebraic decay of g\ (r) (linear 
in the double-logarithmic scale of Fig. [5]). By further 
increasing U, gi(r) still displays a power-law decay but 
falls off more rapidly as a function of r. This is in analogy 
with the homogeneous case [21( where interactions cause 
quantum fluctuations and loss of coherence. At much 
larger values of the interaction energy this trend drives 
the system in the Tonk- Girardeau regime that cannot 
be described by the mean field theory. A similar analysis 
performed at fixed U shows a monotonic reduction of the 
coherence when increasing the disorder strength. 

To better understand the physics behind the functional 
dependence of gi(r), it is useful to focus on the 2-point 
correlation function expressed in Eq. ((5J). Fig. [5] shows 
G(ro, r) for two cases in the Bose glass and quasiconden- 
sate phases respectively. It is evident that the correlation 
in the superfluid phase, although slightly modulated by 
the underlying disorder profile, displays a smooth decay. 
On the other hand, the Bose glass case is characterized 
by abrupt jumps in the 2-point coherence that separate 
relatively coherent zones. It is immediate to notice that 
these jumps coincide with the nodes of the lowest en- 
ergy excitations as can be checked comparing with the 
corresponding vj^ (r)-modes in Fig. [SJd. The size of the 
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FIG. 7: (color online) Averaged D{E) at fixed A = 0.8 E c for 
various interaction energies. The triangles denote two Bose 
glass cases (U = 0.48 — 0.8 E c ), whereas squares and circles 
mark two superfluid phases (U = 1.12 — 1.44 E c ). 



jumps is related to the amplitude of the excitations, that 
increases for decreasing energy. 



V. DENSITY OF STATES AND LOCALIZATION 

Given that the loss of coherence in the Bose glass phase 
is due to low-lying Bogoliubov excitations, we study their 
properties in detail. In particular we are interested in the 
DOS and the localization properties of the v j^-modes, 
that are the only ones playing a role at zero temperature. 
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The DOS is denned as 



D{E)=Y,S{E-E L 



3)1 



(11) 



where the Ej are the positive Bogoliubov energies. This 
quantity is expected to approach a constant value for 
E —t in the supcrfluid phase, as for phonons in ran- 
dom elastic chains [l6| . Moreover, some theoretical stud- 
ies @, Hl| have argued that the low-energy DOS should 
remain constant also in the Bose glass phase. The results 
of our mean field calculation disagree with this latter pre- 
diction in the mean field limit. As shown in Fig. [7l in 
the case dominated by interaction, D(E) approaches a 
constant value at low energy, while it develops a power- 
law divergence in the Bose glass phase. In fact, by de- 
creasing the interaction strength at fixed disorder, the 
DOS for E — > increases and it starts diverging follow- 
ing a power-law beyond the phase boundary. In addition, 
the slope of the power-law increases monotonically going 
deeper in the insulator phase. 

We can connect this divergence to the behavior of 
g\. Inspection of Eq. ([5J shows that at long dis- 
tances, the main contribution comes from the term 
J2j \ v ±j{ r )\ 2 - Thus, we can rewrite G(r, ro) ~ 
exp[- J\v E ±(r)\ 2 D(E)dE] (see also Ref. [13). Here 
|ve_l(t*)| 2 is defined as the local density of Bogoliubov 
excitations per unit energy, i.e. 



J2j S(E - Ej) 



(12) 



This quantity is expected to diverge as l/E, as it is the 
case both in the homogeneous system and in the case of 
weak links between junctions [231 ] (as also checked nu- 
merically for a disordered potential). Hence, the change 
in <7i(r) can be linked to a change in D(E), that changes 
from constant to a power-law divergence for E — ► 0. 
These low energy excitations change the phase between 
weakly coupled neighbouring islands. The DOS is there- 
fore directly related to the statistics of the strength of 
the weak links investigated in Ref. [l7j . Our diverging 
DOS is in agreement with their analysis. 

The discrepancy between our results and former pre- 
dictions @, [HJ could also be due to a different nature of 
the phase transition in the weakly and strongly interact- 
ing regimes, as also suggested by a recent renormaliza- 
tion group analysis [T^. These results might imply two 
glassy phases characterized by different properties. In 
particular the analysis by Altman et al. has found a so- 
called random- singlet phase characterized by a divergent 
DOS and claims that the phase transition at strong dis- 
order belongs to a different universality class with respect 
to the weak disorder transition [6[. The random-singlet 
phase is specific of systems with particle-hole symmetry, 
as it is the case for the Bogoliubov model studied here. 

We now turn to the localization properties of the Bo- 
goliubov modes. As pointed out in the previous sections, 
in presence of interaction, delocalized low-energy phase 



fluctuations are the main mechanism of reduction of co- 
herence. The localization properties of these excitations 
is still under debate [H, [3| ■ As a measure of localiza- 
tion we choose the inverse participation number (IPN), 
that directly gives an estimate of the spatial extent of 
the wave function. We notice however that IPN can be 
very different from other quantities characterizing the lo- 
calization as, for instance, the exponential decay length 
of the wavefunction tails [24[. The IPN is nevertheless 
the most relevant characterization for our purposes. In- 
deed, a phase change over a long distance can only be 
produced by an excitation whose wave function is signif- 
icantly non-zero at points very far apart in space (even 
if it has rap idly decaying exponential tails outside these 
regions [24|). The IPN is defined as 



1 _ fdr\v ±J (r)\ 4 
h (fdr\v ±j (r)\* 



(13) 



and the corresponding realization-averaged quantity is 
L a (E) = J2j IjS(E-Ej)/D(E). In Fig. [8jthe results are 
shown for fixed disorder and varying interaction strength. 
As it can be noticed, the IPN always shows a power- 
law divergence E~ a for E — > 0, with a increasing when 
going deeper in the superfluid phase. The finite size of 
the simulation limits our analysis for large U, where the 
IPN saturates, as can be seen in Fig. [51a. We draw a 
horizontal line to separate the region of the plot that is 
not affected by the finite size of the numerical sample. 

As it is shown in Fig. [5Jb, we find an exponent varying 
continuously around 1 and crossing that value in corre- 
spondence with the boundary computed by the correla- 
tion length. In the quasicondensed phase a > 1, whereas 
we find a < 1 in the Bose glass phase. It decreases 
by lowering the interaction strength, apparently linearly 
vanishing for U — > as predicted for the non interacting 
case, where the localization goes to a constant value at 
low energy. We can compare our results with a recent 
theoretical work [l3j that predicts an exponent a = 1 
at the phase transition and 1 < a < 2 in the superfluid 
phase. We find a full agreement with this prediction, al- 
though it is not possible to reach the case with exponent 
a = 2 because of the finite size of our simulations. 

A precise characterization of a phase boundary via the 
DOS turns out to be difficult, whereas the IPN analysis 
gives the boundary with high accuracy as evident in Fig. 
lb. 



VI. SUPERFLUID FRACTION 

An alternative way to characterize the thermodynamic 
phase of the gas is the computation of the superfluid frac- 
tion. The usual approach to superfluidity is the two-fluid 
picture and the distinction between normal fluid and su- 
perfluid resides in the different response to a small veloc- 
ity field. This turns out to be equivalent to imposing a 
phase twist in the boundary conditions [25| . The super- 
fluid fraction can be computed from the energy difference 
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FIG. 8: (color online) a. Averaged IPN for fixed A = 0.8 E c 
varying U . The cases with U — 0.48 — 0.8 E c are in the Bose 
glass phase, whereas for U — 1.12 — 1AAE C are superfluid. A 
saturation of the IPN occurs at large lengths because of the 
finte size of the simulated system, the horizontal black line 
represents a confidence limit, b. Computed exponents of the 
power-law divergence of L a (E) for different values of U. 

between the system at rest and the moving one, namely 
2mL 2 Eq — E Q 

where G is the total phase twist. Eq is the energy of 
a condensate with twisted boundary conditions {^(L) = 
^(0)e ie ) and E is the ground state energy of the system 
at rest. We consider a small phase twist, (9 = ir/32 <C 
7i"), to avoid excitations and level crossing. With a gauge 
transformation ^f(x) — > ^e l&x ^ L , the twisted boundary 
problem is mapped on a problem with periodic bound- 
ary conditions, with shifted momentum p — > p + HQ/ L, 
so that V — > V + iQ/L. This substitution enters both 
the Gross-Pitaevskii equation and the Bogoliubov-de 
Gennes equations (|3]). In the homogeneous case, Galilean 
invariance ensures that these latter give no contribution 
to the energy difference (26| . whereas in the disordered 
case they have shown to develop a finite contribution, 
as emerged from numerical simulations. In Fig. [9] the 
superfluid fraction is reported as a function of the inter- 



action energy for three different fixed values of disorder. 
The computation of the twisted problem being quite de- 
manding in terms of computational resources, thus lim- 
iting the number of disorder realizations, results in the 
large error bars in the plots. The effect of the finite size 
simulation is that the computed superfluid fraction is ex- 
pected to be slightly larger than the real one and only 
a careful size-scaling analysis gives reliable quantitative 
information about the thermodynamic limit. To improve 
the averaging procedure we compute the energy of a sys- 
tem with twisted boundary conditions splitting it into TV 
bins as 

N 

s e = X>*(0<-0i-i) a , (15) 

i 

where psi is the superfluid fraction of the i th bin and 
9i,0i-i are the phases at the boundary of the i th cell. 
With the constraint 

£> = e, (i6) 

i 

it can be shown that minimizing the energy (|15p corre- 
sponds to taking the harmonic average of the superfluid 
fractions, i.e 

»- (?£)"'■ 

For this reason each point shown in Fig. [9] is computed 
as a harmonic mean of the superfluid fraction of each 
realization and the error bars are computed accordingly. 
The shaded zone shows the phase boundary predicted by 
studying the long-range decay of the one body density 
matrix. 

In Fig. it is evident that the boundary computed by 
means of the correlation length coincides with the predic- 
tion based on superfluidity. In fact, the zero- value for the 
superfluid fraction is consistent with the error bars of all 
the cases belonging to the insulator phase, whereas the 
superfluid points acquire a finite ps- It is worth notic- 
ing that the average procedure is most demanding when 
close to the boundary and is reflected in larger error bars 
in the proximity of the phase transition. 

VII. PHASE DIAGRAM 

With the methods explained so far we can characterize 
the phase of the Bose gas. We are able to draw the mean 
field phase diagram of the ID Bose gas at zero temper- 
ature as a function of disorder and interaction energies. 

We presented an earlier version of this phase diagram 
in Ref. [l||. In Fig. [TO] we show an extended phase 
diagram that includes the TF regime. It clearly shows 
two different trends depending on the ratio between the 
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FIG. 9: (color online) Superfluid fraction for 3 values of disorder: a. A = 12.8 E c , b. A = 0.8 E c , c. A = 0.016 E c . The 
average superfluid fractions and their error bars are shown as a function of the interaction energy. The shaded zones mark the 
phase transition computed through the degree of coherence. 
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FIG. 10: (color online) Sketch of the phase diagram of the ID 
Bose g £LS EIS 3, function of interaction and disorder. (A): Bose 
glass; (O) 1 quasi-condensate. 



characteristic energies at the transition and E c . These 
regimes can be identified by the ratio k = U / E c , in fact 
the limit k«1 represents the WN limit, where the heal- 
ing length is much longer than the disorder correlation 
length. The opposite case k ^> 1 marks the TF regime. 
The numerical results give two power-law dependencies of 
the boundary in these limiting cases A/E c = C(U / 'Ec) 1 , 
with 7 equal to 3/4 and 1 respectively. The lower part 
of the phase diagram represents the WN limit: in this 
regime a sing le energy scale characterizes the disordered 
potential [27| (see Appendix) 



1/3 



Thus, assuming that at the transition the interaction en- 
ergy is proportional to Eq, we directly obtain 



E r V E r 



4/3 



(19) 



(18) 



that correctly reproduces the power-law found numeri- 
cally. In the opposite regime, i.e. the TF regime E c , 
the scale of the potential is much larger than the typi- 
cal length involved in the modulation of the ground state 
wavefunction (£ <C rj). In this limit E c is no longer rel- 
evant and the thermodynamic phase is only determined 
by a one-to-one competition between disorder and inter- 
action energies. The computed phase diagram is valid 
in the weakly interacting limit and becomes exact in the 
limit p — > oo at constant gp. The phase diagram shows 
an infinite slope in the origin, in fact the power-law is 
smaller than 1, and this is in agreement with previous 
theoretical calculations [ll[. This condition implies that 
in an experiment where the correlation length r] is re- 
duced at constant interaction and disorder amplitudes, 
one would always end up in the superfluid phase. The 
obtained phase diagram is strikingly similar to the dia- 
gram obtained by investigating the change in the density 
profile of the II? Bose gas [28|]. A precise link between 
the phase transition and the fragmentation of the ground 
state wavefunction will be the object of a future study. 

For the proportionality constant in the WN limit, we 
numerically find C\ ~ 1.1. This is in good agreement 
with the prediction [29j that the proportionality between 
Eq and U should be approximately 1, (Eq/U ~ 1). 

As stated above the mean-field description does not 
hold in the strongly interacting regime, that leads finally 
to an interaction dominated Bose glass phase. Conse- 
quently, the reentrant Bose glass phase obtained in the 
discrete model for strong interaction @ cannot be de- 
scribed within this mean-field model. Moreover, the 
model cannot be applied in the disorder-dominated case, 
where the coherence extends only within a few maxima 
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of the density and the system is in the so-called Lifshitz 
glass phase 
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VIII. CONCLUSIONS 

We have studied the phase diagram of a ID Bose-gas 
at zero temperature in presence of correlated disorder. 
We analyzed the changes in the Bogoliubov excitations 
that entail the phase transition: we have found that the 
DOS diverges in the Bose-Glass phase while it approaches 
a constant value in the quasicondensed case. Moreover 
the localization of the excitations always shows an E~ a 
divergence and a = 1 marks the phase transition. We 
have established the phase diagram by inspecting the 
long range decay of the one-body density matrix. This 
analysis led to the identification of two regimes in which 
the boundary follows a power-law relation between disor- 
der and interaction: a WN zone, where a 3/4 power-law 
relation holds, and a TF regime, where the relation be- 
comes linear. This phase diagram has been confirmed by 
inspecting the superfiuid fraction of the system. 

This theoretical analysis could be very useful for future 
investigations aimed at the determination of the super- 
fluid to Bose glass phase transition in ID weakly inter- 
acting alcali gases. 
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APPENDIX A: WHITE NOISE LIMIT 

The WN limit argument applies to arbitrary dimen- 
sionality D to predict the behaviour of the phase bound- 
ary at the origin. The correlation of the potential is given 
by [30( (V(r)V(r')) = A 2 / r _ r /, where for a gaussian cor- 
relation we have / r - r ' = e~( r ~ r ) l 2r) . In the WN case 
the potential is equivalent to a delta correlated one, 



(W(r)W(r')) =wS(t-t'), 



(Al) 



so that, if we assume h 2 /2M = 1, w has the dimen- 
sionality \E 2 ~ D I 2 \ and in the WN limit every quantity 
having dimension of an energy must be proportional to 
Eo = w 2 ^ i ~ D ' 1 . On the other hand, if we take the 
limit for r\ — > of the Gauss-correlated potential we get 
/ r _ r ' oc rj D 5(r — r'). Comparison with Eq. (|A1|1 implies 
w ~ A 2 r) D . In particular the critical interaction energy 
should be proportional to E in the WN limit. Using 
rf cx E c 2 , we conclude that 



(A2) 



U 

E c \E C 

The prediction for the 2D and 3D cases are respectively 
K = 1/2, 1/4. Thus the slope in the origin remains infi- 
nite and the difference with respect to the linear relation 
in the TF limit becomes more pronounced in higher di- 
mensions. These results are the same that have been 
found by Falco and coworkers fill ]. 



[1] P. W. Anderson, Phys. Rev. 109, 1492 (1958). 

[2] J. Billy, V. Josse, Z. Zuo, A. Bernard, B. Hambrecht, 

P. Lugan, D. Clement, L. Sanchez-Palencia, P. Bouyer, 

and A. Aspect, Nature 453, 891 (2008). 
[3] G. Roati, C. D'Errico, L. Fallani, M. Fattori, C. Fort, 

M. Zaccanti, G. Modugno, M. Modugno, and M. Ingus- 

cio, Nature 453, 895 (2008). 
[4] B. Deissler, M. Zaccanti, G. Roati, C. D'Errico, M. Fat- 
tori, M. Modugno, G. Modugno, and M. Inguscio, 

arXiv:0910.5062vl. 
[5] Y. P. Chen, J. Hitchcock, D. Dries, M. Junker, 

C. Welford, and R. G. Hulet, Phys. Rev. A 77, 033632 

(2008). 

[6] T. Giamarchi and H. J. Schulz, Phys. Rev. B 37, 325 
(1988). 

[7] M. P. A. Fisher, P. B. Weichman, G. Grinstein, and D. S. 
Fisher, Phys. Rev. B 40, 546 (1989). 

[8] N. V. Prokof'ev and B. V. Svistunov, Phys. Rev. Lett. 
80, 4355 (1998). 

[9] S. Rapsch, U. Schollwock, and W. Zwerger, Europhys. 
Lett. 46, 559 (1999). 
[10] L. Fallani, J. E. Lye, V. Guarrera, C. Fort, and M. In- 
guscio, Phys. Rev. Lett. 98, 130404 (2007). 



[11] G. M. Falco, T. Nattermann, and V. L. Pokrovsky, Phys. 

Rev. B 80, 104515 (2009). 
[12] P. Lugan, D. Clement, P. Bouyer, A. Aspect, M. Lewen- 

stein, and L. Sanchez-Palencia, Phys. Rev. Lett. 98, 

170403 (2007). 

[13] V. Gurarie, G. Refael, and J. T. Chalker, Phys. Rev. 
Lett. 101, 170407 (2008). 

[14] N. Bilas and N. Pavloff, Eur. Phys. J. D 40, 387 (2006). 

[15] L. Fontanesi, M. Wouters, and V. Savona, Phys. Rev. 
Lett. 103, 030403 (2009). 

[16] T. A. L. Ziman, Phys. Rev. Lett. 49, 337 (1982). 

[17] E. Altman, Y. Kafri, A. Polkovnikov, and G. Refael, 
arXiv:0909.4096vl. 

[18] C. Mora and Y. Castin, Phys. Rev. A 67, 053615 (2003). 

[19] S. Shevchenko, Sov. J. Low Temp. Phys. 18, 223 (1992). 

[20] V. N. Popov, JETP Lett. 31, 526 (1980). 

[21] L. Pitaevskii and S. Stringari, Bose-Einstein Condensa- 
tion (Oxford University Press, 2003). 

[22] P. Nisamaneephong, L. Zhang, and M. Ma, Phys. Rev. 
Lett. 71, 3830 (1993). 

[23] G. S. Paraoanu, S. Kohler, F. Sols, and A. J. Leggett, J. 
Phys. B 34, 4689 (2001). 

[24] B. Kramer and A. MacKinnon, Rep. Prog. Phys. 56, 1469 



10 



(1993). 

[25] E. H. Lieb, R. Seiringer, and J. Yngvason, Phys. Rev. B 

66, 134529 (2002). 
[26] I. Carusotto and Y. Castin, Comptes Rendus Physique 

5 (2004). 

[27] B. I. Halperin and M. Lax, Phys. Rev. 148, 722 (1966). 



[28] P. Lugan, PhD thesis. 

[29] I. Aleiner, B. Altshuler, and G. Shlyapnikov, 

arXiv:0910.4534vl. 
[30] V. Savona and R. Zimmermann, Phys. Rev. B 60, 4928 

(1999). 



